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Abstract 

Correspondence in quantum chaotic systems is lost in short time scales. Introducing some noise 
we study the spectrum of the resulting coarse grained propagaor of density matrices. Some differen 
methods to compute the spectrum are reviewed. Moreover, the relationship between the eigenvalues 
of the coarse-grained superoperator and the classical Ruelle-Pollicott resonances is remarked. As 
a concequence, classical decay rates in quantum time dependent quantities appear. 
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I. INTRODUCTION 



Since the beginings of quantum theory, the problem of the transition from the quan- 
tum to the classical world has posed interesting and challenging questions. The concept of 
decoherence []| has introduced some clarity, yet many questions remain to be answered. 

Quantum chaos studies the emergence of classical chaotic behavior in quantum systems as 
h — > 0. The introduction of Gutzwiller's trace formula|2| and the developement of periodic 
orbit theories and in addition the understanding of the influence of random matrix 

theories^ gave quantum chaos a boost of interest. One key feature of quantum chaotic 
systems is that correspondence is lost in very short time scales (Eherenfest time (?]]). The 
introduction of noise (in a controlled way) restores correspondence and allows to study it 
from the point of view of the appearance of classical properties from the quantum systems. In 
this way, the classical Lyapunov exponent plays an important role in the medium time decay 



3,0] and the Loschmidt echoQ, U| 



of time dependent quantitites like the linear entropy 

The study of this correspondence for long times, like for example the asymptotic behav- 
ior of correlation functions, requires the understanding of the spectral properties of both 
quantum and classical propagators and their relationship in the h — > limit. 

In this review the latter approach is discussed in detail although the former is also ad- 
dressed. Some recent techniques to compute the spectrum of both the classical and the 
quantum propagator of densities are described. We describe in detail some efficient methods 
used by the authors recently^, 12, Q. The systems are quantum maps on the torus. By 
introducing diffusion in the form of Gaussian noise we obtain a coarse grained quantum 
propagator |3q. The Gaussian noise introduces a natural truncation that is understood us- 
ing phase space representation theory (chord funtion and Wigner function). In addition, 
we develop an iteration method that efficiently provides the leading part of the spectrum. 
Relation to the classical Ruelle-Pollicott resonances is then established by taking the limit 
h — > and also by studying the long time decay rates of several quantities. 

This review is organized as follows. In Sec.|TI]we define the phase space we will work with 
and summarize some operator theoretic properties of representations in phase space. A brief 
introduction to the theory of Open Quantum systems is provided in Sec. IIHI In sections IIVI 
and |3 quantum maps and the model used for open quantum maps are described. The kind 
of Gaussian noise used, in term of phase space translations and how it acts when composed 
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with a unitary superoperator, yielding a coarse grained propagator, is developed in Sec. IVII 
while the issue of correspondence, specifically through the relations between classical and 
quantum spectrum is addressed in Sec. IV111 Finally in Sec. I Villi some numerical methods 
to compute the spectrum are described. The consequences for the asymptotic decay for long 
times are discussed in Sec. IV111DI 



II. PHASE SPACE 



A. Torus quantization 

All the calculations are made for systems with a bounded phase space, which for simplicity 
we take to be a 2- Torus (T 2 ). 

There are some well-known features that appear in the quantization of T 2 . Since the area 
is finite, and normalized to unity, the Hilbert H,n space is finite and of dimension N. The 
value of N determines an effective Planck constant 

h=— . (l) 

The position and momentum bases are discrete sets 

B q = {\ qi ), qi = 0,l,...N-l} 
B p = {\p J ),p J =Q 1 l ) ...N -1} 

which define aniVxiV grid G^. We identify with |<&) the position qi/N and with \pj) the 
momentum Pj/N. They are related by the discrete Fourier transform (DFT) of dimension 

N, 

(p\q) = -f=e-W^>. (3) 



Pure states are represented by iV-dimensional vectors |0) in Tijq. Mixed states are best 
described by the density operator p. If we identify Liouville space TLn 2 with the space of 
comlplex N x N matrices, then density operators associated to states mH,N form a subset 
of TLn 2 of self-adjoint, positive semidefinite matrices, with unit trace. Moreover 

Tr(p 2 ) < 1, (4) 

with the equality holding iff the state is pure. 
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We now define the translation operators on T 2 . In M 2 the Weyl group defines the trans- 
lation operator of q\ in position and p\ in momentum, as 

f ig>p) = e -m(nP-PiQ) . ( 5 ) 

Translations on the torus are somewhat harder to define because infinitesimal operators 
Q and P with the canonical commutation rules cannot be defined in a discrete Hilbert space. 
Nevertheless, finite cyclic shifts U and V such that, if \q) and \p) are position and momentum 
eigenstates respectively, can be defined [lJ] as 

U<n\q) = \q + qi ); U^\p) = e -(^/N)gip\ p ) ; 

(6) 

VK\p) = \p+pi); V pi \q) = e (2 " i/7V)pi9 |g) 
These operators satisfy the comutation rule 

Symmetrizing Eq. (J7J) one gets the transaltion operators on the quantized torus 1 as a com- 
position of position and momentum shifts, 

= yp(j<i e -{™/N)qp ^ 

with q,p = 0, . . . , N — 1. There are N 2 operators T^ q ^ and have the following property 

Tr(flf p ) = N5 a p, (9) 

for all a,f3 (when convenient we use greek letters to rperesent a phase space point, e. g. 
a = (q,p)). Therefore, they constitute a complete orthogonal set in the Hilbert-Schmidt 
inner product. Moreover, they have the group composition rule 

T a Tp = e^ N ^f a+p (10) 

where a A (3 is the usual wedge product. 



1 We do not distinguish the notation of translations on M 2 and T 2 because we mainly use the latter. 
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B. Phase space representations of states and operators 



The formulation of quantum mechanics in phase space has proven useful in the under- 
standing of correspondence because phase-space is the natural stage for classical mechanics. 
If one looks for a quantum analog of a probability distribution in phase space the Wigner 
function appears to be the best choice. If p is a density operator in a continuous Hilbert 
space the Wigner function is usually defined as 

/°° dr 
J*/*{ q - x /2\fi\ q + x /2). (11) 
■oo 2™ 

Although it is not necesarily positive everywhere, this function is real valued and the integral 
along any line in phase space yields the correct marginal distributions. 2 Bertrand and 
Bertrand showed that this properties toghether with 

Tr(p!p 2 ) = 2nh J W X W 2 da, (12) 

uniquely determine the Wigner function. 

Quantum discrete phase space descriptions can also be made. There are many works on 
different subjects that take advantage of this line of research. 3 We briefly review general 
operator theoretic notions of representations. 

Operators in T~Cn 2 can be represented by their c-number function in some operator basis. 
There are N 2 linearly independent operators in Ti^ with the Hilbert-Schmidt inner product. 
Let {S a }a=o 1 be a complete set of operators such that 

Tr(St^) = 5 a0 (13) 

then an arbitrary operator B can be expanded as 

N 2 -l 

B=J2 M<*)Sa- (14) 

«=0 

The coeficients of the expansion are 

S B (a) = Tr(SlB) (15) 



2 Measuring along different lines is the basis of a widespread technique such as quantum state (and process) 

tomngraphv 



See 



iiliali3lialiai2o|, to list just a few. 
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and Sb is called the symbol of B in the S a basis representation. 

If the basis operators are the translations T a then we have the characteristic function 
representation 

a 

with the symbol 

C B (a)=Ti(flB), (17) 

The function Cb{ol) is the so-called chord function\ll\ of B, and it owes the name to its 
particular geometrical features. 

The translations defined in Eq. (JSJ) are not periodic in the N x N grid, although they 
form an orthonormal set. If we extend the grid GV to a 2N x 2N grid G^at then the DFT 
of the translations in this grid 

(2A0 2 -1 

A ? = j2Ny E te- (27ri/2JV) ^ (18) 

defines the discrete phase-space point operators. There are 4N 2 of these operators. However 
it can be shown [19] that only N 2 are linearly independent. The expansion 

B = N WB(a)A a , (19) 

a€G 2 N 

determines the discrete Weyl (or center) representation of B. The symbol W p of a density 
operator p in the Weyl representation is the discrete Wigner function (DWF). Miquel, et al. 
[19j show that the DWF defined by Eq. (jTHj) has all the desired properties analogous to those 
of the continuous Wigner function. 

III. OPEN QUANTUM DYNAMICS 

When studying realistic models of quantum systems one usually thinks in terms of the 
specific system of interest (labeled s) and a reservoir (labeled r) with which the system in- 
teracts. Their evolution as a whole is unitary and (we assume) determined by a Hamiltonian 
like 

H = H s + H r + H sr (20) 
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where H s and H r describe the evolution of the system and the reservoir independently and 
H sr is the interaction Hamiltonian. The state of the system is given by the density operator 

f) sr G Ji s ® H r . 

Standard treatment of quantum systems in the presence of noise, generally in the form of 
a reservoir or bath of harmonic oscillators, involves tracing out from the unitary evolution 
of the whole (system plus environment) only the reservoir degrees of freedon yielding a 
nonlinear partial differential equation or master equation. When the bath is assumed to 
fulfill de Born-Markov approximation, which requires that the relaxation times of the bath 
are much shorter that any characteristic time of the evolution of the system, the master 
equation can be put in the form 

% = - l j\H,p}+ l -Y, ca { [iU, F]] + [4 pF}] } , (21) 

u 

where p s is the density matrix of the system and Fj are system operators. When the matrix 
of coeficients Qj (GKS matrix (2^) can be diagonalized Eq. (f2Tj) takes a simpler, so-called 
Lindblad 22], form. Gorini, et al. [21] show that equations of the form Eq. (|2~T]) generate 
a one parameter dynamical semigroup of trace-preserving, completely positive quantum 
operations. Complete positivity means that a trivial extension of the quantum operation to 
a larger Hilbert space yields a positive operation. 4 

A trace preserving, completely positive quantum operation S can always be written in a 
generalized operator sum or Kraus [2j| form 

S(p)=J2^W,pMl (22) 

where c(/i) > and are system operators. If N is the dimension of the state space, then 
p is no bigger than iV 2 . The map can be shown to be TP if 

Y,c{p)MlM, = i. (23) 

IV. QUANTIZED MAPS 

We will focus our attention on discrete time systems, i.e. maps. Instead of modelling 
interaction with the environment by proposing the system operators Pi that appear in the 

4 See Preskillji^ or Nielsen and ChuangQ 
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master equation, we directly model the effect the noise has on the system by giving the noise 
superoperator specifically in Kraus operator sum form. 

General dynamical are defined by a first order differential equation 

x = /(x) (24) 

where x is a point in phase space. Usually the solution (and the visualisation) is simplified 
by means of a discrete map which arise by observation of the map at fixed time intervals 
(stroboscopic map) or when the flow passes transversely through some hypersurface V in 
phase space (Poincare section). So now time is discrete, and in general labeled n, and the 
equation is a recurrence, 

y n+1 = M(y n ), y G V. (25) 

Area preserving maps are characterized by a finite canonical transformation. The uni- 
tary representation of this canonical transformation [26] is what is usually considered the 
quantum version of the map. Although there are not general procedures to obtain this rep- 
resentation, in the real two-dimensional plane R 2 , the procedure is relatively standard. |27j] 
An aproximation of the unitary propagator in the semiclassical limit can be written as 

% 



i d 2 S ^ 1/2 



U ^ )=[ hd^ 2 ] 6XP 



h S(q 2 ,qi] 



(26) 



where S is the action along the unique classical path from q\ to q 2 , and where for simplicity 
we do not consider the existence of multiple branches and Maslov indices. Only for linear 
symplectic maps on R 2 this unitary propagator is exact, and then S is minus the quadratic 
generating function of the linear transformation. By different ad-hoc procedures, some of the 
best known classical maps have been 

and shears, and also chaotic, such as cat maps, Q baker maps,j3, Q and the standard 
map. |3l| Moreover, noncommuting nonlinear shears, or "kicked" maps, can be quantized as 
well as time dependent periodic Hamiltonians. 3^] 

On T 2 a quantized map is a unitary operator U G Ti^ 2 - The evolution of density matrices 
by this map is given by 

p n+1 = Up n Ul (27) 

The last equation can be written as a unitary map U : H-n 2 ~~ * ^tv 2 i n Liouville space, that 
takes density matrices to density matrices 

pn+i = (UQuAp n 
8 



= ua 



(28) 



We introduce the symbol to represent adjoint or left-right action, 

((J U^j p = Ad(U)p = Uptf (29) 

Maps of density matrices to density matrices are called superoperators or quantum opera- 
tions. 5 Further properties of more general (i.e. non-unitary) superoperators are given in the 
subsequent sections. 



V. OPEN QUANTUM MAPS 

We propose as a simplification to model discrete time open quantum evolution with 
a nonunitary propagator (superoperator) S that consists of two steps. First the unitary 
U step corresponding to noiseless evolution. Then the noise acts through a nonunitary 
superoperator D e , where e characterizes the coupling strength. Thus the evolution given by 
the two-step noisy propagator looks like 

p n+1 = S e (p n ) = D e o U(p„) (30) 



This type of two-stage schemes have been studied recently in several works. |9lll2t 113. 



30, [3 



33 



34. 



351 ] They appear naturally in situations in which the noise is negligible during the unitary 
evolution. Suppose the solution comes from the direct exponentiation of the Lindblad equa- 
tion. Then if to leading order, the unitary part and the noise part commute, then the two 
parts can be separated into a product of two propagators. One clear physical example is the 
micromaser 



36( . There are some other cases like a kicked map in which the interaction takes 
place between kicks. The case of a billiard inside a bath, if the interaction with the walls 
is cosidered essentially unitary while free propagation is noisy, constitutes another example. 
Another example could be an iterated quantum algorithm, if the errors due to imperfect 
gates can be neglected but the channel is noisy. 



We use indistinctly both denominations. 
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VI. DIFFUSIVE NOISE AND COARSE-GRAINED QUANTUM PROPAGATOR 



We now proceed to describe the type of noise D considered for the two step propagator 
of equation Eq. (J30j) . We will use noise models given explicitly in Kraus form 



i 

with U a unitary operator and pi > 0. From Eq. ()23|) we need 

5> = 1 



(31) 



(32) 



for trace preservation. This condition also implies that the sum in Eq. 1)31)1 is a convex sum 
of unitary operators. Any map that can be written as a convex sum of unitary operators 
is a random unitary process (RUP). The physical interpretation of a RUP is simple: the 
unitary operation UipUj is applied to the state p with probability pi, and the sum inEq. ()31j) 
provides a simple model for a quantum Langevin-type process. It is easy to show that a 
trace preserving RUP is also unital or identity preserving, that is 



D(J) = I 



which means that the eigenvalue 1 has the maximally mixed state 



Poo = I/N 



(33) 



(34) 



as eigenstate. Maps that are both TP and unital are called bi-stochastic and they are 
contractive maps which means that the spectrum contained in the unit circle. 

We propose as diffusive noise the following RUP which is diagonal in the chord represen- 
tation 

D « = ^£<U°) f « £! ( 35 ) 

a 

with C e (a) a Gaussian of width proportional to e. The chord function of D e is 



Tr[fjD e f;] = Tr 



a 



Tr 
1 



(36) 
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where C e (3) is the DFT of C e (a) and is also a Gaussian of complementary width oc 1/e. As 
a consequence if the chord function of a state p is say g(a), the action of D e in the chrod 
representation is 

g'((3)=C e ((3)g((3) (37) 

and the long chords {(3 3> 1/e) are supressed. The chord function is related to the DWF 
by a Fourier transform so the long chords are related to the small structures in the wigner 
function. The small scale structures in the DWF are associated to quantum interference 
and coherence. Therefore erasing the long chords is equivalent to erasing short wavelength 
oscilations in the DWF resulting in a decoherent process. 

In analogy with Eq. (J37|) the chord re [presentation provides a picture of what the effect of 
the noise when composed with a unitary process. Composition with a unitary map U Xft 
yields again a RUP (and thus a contractive map). Let S = D e oU, with U = U (ft. Then 
in the chord representation it is just 

S(aP) =C e (a)U(a[3) (38) 

where U(a/3) = Tr[T^UT^(ft], and Therefore the Gaussian function C e (a) modulates the 
chord funtion of U supressing the long chord components. This is in fact a coarse graining 
of U and we shall call S e coarse-grained quantum propagator. 



VII. CHAOS AND CORRESPONDENCE 



Let us first give the now standard definition of quantum chaos. By quantum chaotic 
systems we mean those whose classical counterpart is chaotic. Usual approaches to under- 
standing chaos and quantum mechanics range from the concept of universality in random 
matrix theory I™, and also since the celebrated Gutzwiller|2l trace formula the use of periodic 

nn 

orbits in the computation of the eigenvalues and eigenf unctions |4J, |5( became fundamental. 

Recently yet another different approach has been explored. For strongly chaotic (Axiom- 
A 6 ) systems correlation functions, as well as any time dependent quantity decay exponen- 



tially. The decay rates are closely related to the spectrum of t 



phase space. The idea is to establish, be it numerically 



34. 



re propagator of densities 



, experimentally 
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6 Axiom-A or Anosov-A are uniformly hyperbolic systems with strong mixing. 
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or theoretically[35f, a relationship between the classical decay rates and the spectrum of the 
noisy quantum propagator. 



Ruelle-Pollicot Resonances 



Let M be an invertible area preserving map. 

Formally the evolution of a distribution p(x) in phase space is given by the Liouville 
equation 



dp 
dt 



Cp = {H, p} 



(39) 



where the Liouville operator £ is the generator of the Perron Frobenius (PF) operator V 



(40) 



When possible we can reduce the dynamics to a stroboscopic area preserving map M and 
the solution in terms of the PF and M can be ritten asjjj 



Pn+l = Vp r , 



dx'5(x-M(x'))p„(x'). 



(41) 



The PF is unitary in the space of square integrable functions L 2 . 

Effective irreversibility in chaotic systems manifests in the exponential decays of corre- 
lation functions. For two functions /, g in phase space we will use the following simplified 
expresion for the correlation function 



C t (f,g) = (f,V t g)-(I,f)(I,g) 



(42) 



where ( , ) denotes inner product and I is the uniform density. For strongly mixing systems 
Ct{f,g) goes to zero as t — » oo. The decay rates are the so-called Ruelle-Pollicott(RP) 
resonances. For Anosov A systems they are the poles of the resolvent of the PF in a higher 



Riemannian sheet 



39|. There are different methods to unveil the resonances from the PF. 



One method is directly interpolating the exponential decay in the time evolution 40]. Other 



more sophisticated consist in Fredholm determinant techniques 



411 ] and there are also cycle 



expansion [42|] and trace formula approaches [3 
study of decay of correlations by Brini et al. 



Markov Partitions were introduced in the 



The RP resonances can also be related to the spectrum of V acting on a functional space 
different from L 2 . Recently Blank et al. [44] provided an accurate theoretical description 
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lm{X) 




Re(X) 



Figure 1: Schematic representation of the spectrum of V . The dots represent the point spectrum 
(<T ps ) while the circle of radius r represents the essential spectrum. 

of the spectrum of V for Anosov A systems. The idea is that for area preserving maps 
if V is restricted to spaces like L p where the norm essencially provides information about 
size, the spectrum does not provide much information about the dynamics. Therefore they 
use strongly anisotropic Banach spaces of functions especially adapted to the dynamics 7 
composed of asymmetrical sets of functions which are smooth along the unstable direction 
and increasingly singular along the stable direction. They showed the operator restricted 
to this space to be quasi- compact. Moreover they showed that the spectrum is stable under 
pertubation and proposed a discretization scheme to compute the spectrum from finite size 
matrices. The spectrum thus found (represented sechematically on Fig. ^) can be described 
as follows: a point spectrum, with finite multiplicity 



for some radius r > 0, and an essential spectrum inside the circle of radius r. Notice that 
for hyperbolic systems the gap between Ao = 1 and Ai is strictly non-zero. This is essential 
for exponential decays to occur. Now assuming that the eigenvalues are nondegenerate (and 
that r < 1) then we can make an approximate spectral decomposition of V as 



ps — 



{Ai,r < A, < A = 1} 



(43) 




(44) 



7 Similar spaces were also used by Rugh[4i 



13 



where 7Zi, are the right (left) eigenfunctions. Using this spectral decomposition explic- 
itly in Eq. (J4*2~j) we get an expansion of Cf g like 



Cfg{t) 



(45) 



clearly dominated by the point spectrum <7 ps . From Eq. ()45j) it is clear that in the limit 
t — > oo we have Cf g (t) ~ A* (/, TZ\)(Ci, g), and the importance of the largest (in modulus) 
RP resonance becomes evident. The preceding analysis allows us to identify the complex 
eigenvalues Aj in the point spectrum of a sp of V with the RP resonances. 

There are many methods (in use lately) to compute the classical PF spectrum. We 
briefly describe two of the most popular. The first one is a truncation met 



example by Manderfeldei al. 



Khodas et al. 



46], 



rod used for 
Hasegawa and Saphir 4?| and was 



nicely summerized by Fishmanj48| in the following way: 



1. Introduce a 
polynomials 



?asis ordered by increased resolution. Typically Fourier modes, Legendre 
47 ] or spherical harmonics 



2. Truncate to a desired dimension N leaving the rapidly oscillating components out. 
This acts as coarse graining of V. 

3. Compute the eigenvalues of the truncated operator and as N is increased observe that 
there is a part that freezes and is independent of N. 

4. Compare the frozen eigenvalues with the resonances obtained by analytic continuation 
of the resolvent. 

It should be duly noted that this method was successfully implemented also for systems with 



a mixed phase space 



49] 



The second method consists of implementing a coarse graining 
it with a Gaussian kernel of width e 



dpf 

V f = V f o v. 



35 



of V by convoluting 



(46) 



Essentially this is achieved by just replacing the delta functions in Eq. (|41|) by normalized 
(and periodized in the case of the torus) Gaussians of width e. This kind of coarse graining 
is equivalent to adding diffusive noise in the equations of motion, [^l] In order to perform the 
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calculations a discretization of phase space is needed. The RP resonances are the eigenvalues 
of the coarse grained PF, V e , after doing an extrapolation of N — > oo and e — > 0. 8 



B. Spectrum of the Noisy Quantum Propagator 



In a recent work Nonnenmacher 
simplified form): 



351 ] proved the following theorem (which we state in a 



Theorem ( Nonnenmacher 35 1). For any smooth map on the torus and any fixed coarse 
graining parameter e > 0, the spectrum of the coarse grained propagator S e = D e o U, 
converges in the classical limit N — > oo to the point spectrum of the classical coarse-grained 
propagator V 6 = T> e o V . 

Where both coarse grained propagators are like the ones described in the previos section 
and the noise is assumed to be Gaussian with characteristic width e. 

Thus the introduction of noise in the quantum propagator makes the classical decay rates 
(RP resonances) emerge naturally. The problem which persists is the need to explore the 
large N region because then the propagator of density matrices is a, much larger, N 2 x N 2 
matrix. However the decays will be ruled by the largest eigenvalues (in modulus) so in fact 
only a small group of the leading of the eigenvalues is really relevant. In what follows we 
describe some methods used recently to compute the relevant part of the spectrum. 



VIII. NUMERICS 



The goal is to explore regions of N "large". Thus the matrix to diagonalize, in principle, 
is N 2 x N 2 . In the large N limit this may represent a problem in terms of time and 
memory requirements. However, as was previously pointed out, only the largest eigenvalues 
(in modulus) play a significant role in the long time bahavior of observables. Therefore 
approximation methods can be used to significantly reduce the eigenvalue problem. 

In what follows we describe three different methods to compute the leading spectrum. 
The first two are described briefly. They have the advantage that the only limitation on the 

8 The limits are non-commutative, i.e. if e — > faster than TV — > oo then the eigenvalues do not "freeze" 
and the spectrum becomes unitary again. However no general, and optimal, relation e(N) has been 
found. 
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number of eigenvalues to compute is the size of the matrices to digonalize. The third method, 
described in more detail, has the advantage that the first few (typically 10) eigenvalues can 
be computed with very small matrices (typically 10 x 10). 



A. Truncation Method: 



Again, the truncation method described in Sec. lVII XI was succesfully used by Manderfeld 
et al. to obtain a nonunitary, coarse grained quantum propagator and to show that classical 
and quantum evolutions look alike (in the N — > oo limit) and that decay rates can be 
computed from the truncated propagator. However, this approach can have problems of 
interpretation if not of physical realizability. Spina and Saraceno 52j recently showed that 
when the quantum propagator is truncated in the way described, which is essentially a 
sharp truncation, the whole operation, unitary evolution followed by noise is not completely 
posititve and cannont be obtained as a quantum operation. 9 



B. Chord Function Method: 

On the other hand using the diffusive noise of Sec. IV1I and the chord representation a 
"smooth" truncation can be implemented (see Spina and Saraceno [52] and Aolita, et al. 
Q]). Indeed the complete positivity is directly related to the positivity of the coeficients 
C e (a). Now since C e (a) is a (periodic) Gaussian then the spectrum C e (j3) is also a Gaussian 
of complementary width (proportional to 1/e). Therefore a large part of the elements of S e 
in the chord representation will be very small. If one sets for example the elements such that 
(3 is larger than 1/(2ttN) by some factor (bigger than one), exactly equal to zero, then, even 
though there will be oscillations in C e (a) (because it is the inverse DFT) they can be made 
negligibly small just by moving the truncation parameter. Using this method the matrices 
to diagonalize can be reduced to a dimension of order N x N and the precision obtained for 
the resonances is very good. The compromise lies in that using very small matrices to get 
the RP resonances means using very strong noise. So N x N is kind of a lower bound in the 
size of the eigenvalue problem to solve using the chord function method. 

9 There can be some debate on weather complete positivity is an essential requirement. We favour the 
opinion that it is a necesary requirement. 
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C. Iteration Method 



The iteration method is a kind of variation of some known iteration methods 



50, 551 e. 



g. Lanczos iteration method|54|). The main requirement is that the spectrum is of the form 
1 > I Ai| > |A 2 | > ... > |Afc| > .... The convergence depends strongly on the size of the gap 
(1 — | Ai|). Consider the two sets of iterated states 

Bf = {Po,Pi, • • • ,Pk, ■ ■ •}; Bb = {po,P-i, ■ ■ ■ ,/3-fc, • • •} (47) 

where /?o is an arbitrary initial state and 

p k = S k J ; p~ k = S\ k p (48) 

Then we write a kind of projection of S e onto the subspace spanned by this (non-orthogonal) 

sets 

[SJy = Tr^S^-) (49) 
and solve the generalized eigenvalue equation 



Det 



Se]y-A[0] y =0 (50) 



where [0]y = Trfpt^-] is the matrix of overlaps. Now if the spectrum has the required 
form, then both [S e ]y and [0]ij decay very rapidly with i,j so depending on the precision 
required a truncation size k can be chosen. Typicaly in 

aaithe conations .adeBia, 
for N between 50 and 450, the matrices diagonalized where of sizes from 10 x 10 to 15 x 15. 

with this method are displayed. The shaded region 
indicates the region in the prameters iV and e inside which the spectrum stabilizes and can 
be considered classical. Care should be taken when approaching the limits iV — > oo and 
e — > because this limits do not commute. If e goes to zero faster than iV to infinity then 
unitarity is recovered. There seems to be a kind of optimal function e(N) to approach the 
classical limit. The fit is not easy to compute because it is very much system dependent. 

In order to show that the values computed are the desired first k largest eigenvalues we 
assume, for simplicity, that the spectrum is non-degenerate. The propagator can be written 
in its spectral decomposition form as 

S e p = J2 Wr[Ljp] (51) 
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Figure 2: Spectrum of the propagator for the perturbed cat map (k = 0.01) for different values of 
e and N. The dashed line shows approximately the division in the space of the parameters between 
quantum (almost unitary) and noisy (almost classical) behavior. 



where Ri ( Lj) is the right (left) eigenf unction. We can assume also that the overlap between 
the initial state po an d the right and left eigenfunctions is non-zero. Then Eq. (J5U|) can be 
re-written as 

Det[7C T (S e , po, A;)S e /C(S e , p , k) - XK T (S e , p , k)K(S £ , p , k)] = (52) 

where K(S t , po, k) is the Krylov matrix whose columns are the succesive iterates of the initial 
state po U P to order k (which is the truncation size we choose) 

7C(S e , po, k) = [po, S e , po,..., S* _1 p ] (53) 

and T denotes matrix transposition. Then using the expansion of po both in terms of Ri 
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and La it can be shown 121 that the determinant can be written as 



Det[A f SA] = 



(54) 



with 



A 



1 A Aq • • • Aq 1 
1 Ai X\ ■■■ Xt 1 



\ 1 A fc _i A^_ x • • • \\_\ J 



and 



/ 



«oA)(Ao - z) 

«i/?i(Ai — z) 








\ 



(55) 



(56) 



y ••• a fe _ 1 /3 fe _ 1 (A fc _i - z) j 

Now A is a Vandermonde matrix and from m the asumption of non-degeneracy, the determi- 
nant is non-zero. Thus the solution of Eq. (|54|) is just given by the zeros of the determinant 
of S which are exactly the desired first k eigenvalues. Even though we do not discuss preci- 
sion we can say that the efficiency of the method depends mainly on two things. First, there 
has to be a finite gap between 1 and Ai. The second, and very important also, is that the 
map should be easy to implement without constructing the whole N 2 x N 2 superoperator. 



ity distributions such as the Wigner function |5fil. in the 
"\ For isolated systems, the characteristic time for the loss of 



D. Long time behavior 

Exponential stretching (with the rate given by the Lyapunov exponent) is the cause that 
correspondence is lost much faster than for regular systems. This discrepancy is observed 
in phase space quasi-probabi 
form of sub-Planck structure. _ 

correspondence in the Wigner function is usually called Eherenfest |j| time. It is generally a 
very short time and can be expressed as t% = A -1 ln(S/h), where A is the largest Lyapunov 
exponent and S is some characteristic action. 

Zurek and Paz conjectured^] that loss of correspondence in chaotic systems would be 
recovered by introducing decoherence. Moreover they stated that initially the linear entropy 

S n = -ln[Tr[p 2 ]], (57) 
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(which is the natural logarithm of the purity) should grow linearly with n (i.e. time), 
with a slope given by the largest lyapunov exponent of the corresponding classical sys- 
tem. This lyapunov behavior would have to be independent of the coupling strengh, that 
is the strengh of the noise, for a large range of values of the coupling. This statement 
was verified numerically 0, 0]. The lyapunov growth takes place until the Eherenfest time. 
Then the linear entropy saturates asymptotically to the value ln(iV) that corresponds to the 
maximally mixed state p^ — I/N. 

Another quantity of interest with similar short time behavior is the Loschmidt echo. It 
was proposed 0] as a measure of irreversibility and could be used to characterize quantum 
chaos. Many recent works^O, 11] show that the lyapunov regime indeed exists for classically 
chaotic systems. This quantity gained interes recently in the context of quantum information 
because of its close relation to fidelity. 

There is yet another classical regime that can be identified. If we assume that there are 
nondegeneracies then in analogy Eq. (jUj) a spectral decomposition can be made but for the 
quantum noisy propagator 



N' 2 



s e p = ^2 -^ Tr U p 



i=0 



Ri 



(58) 



where Ri and are the right and left eigenoperators of S e , and the eigenvalues are ordered 



by decreasing modulus (starting with Ao 
function behaves asymptotically as 



Then, for example, the autocorrelation 



C n = Tr [p p„] - Tr 



iV 



Pa 



N 2 



A Tr [poPoo] Tr [poopo] + Tr typo 



i=i 



1 

iV2 



A"Tr 



PqRi 



Tr 



L\p] ■ 



(59) 



where Ao = 1 and p^ as defined in Eq. ([34)1 . In Table H] the long time behavior of the 
autocorrelation function, the linear entropy (logarithm of the purity) and the Loschmidt 
echo is summerized. In all cases the logarithm grows linearly with time (n) and the slope 
is determined by Ai. Now taking iV and e in the correct region of Fig. |21 Ai is identified 
with the largest RP resonance of the classical system. Thus another classical property has 
emerged from the quantum system by the addition of noise. 
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Quantity 


Definition 


Aymptotic behavior 


Ln 


Autocorrelation 


Tv[p p n ] - l/N 


~ A™ 


~ n\\ 


Linear Entropy 


-HTr[pl}} 


~ -2nln[Ai] 


2nln[Ai] 


Loschmidt echo 




~ A?(Ai) n 


~ n(Ai + Ai) 



Table I: Long time ehavior of the aurotcorrelation functionC n , the linear entropy S n and the 
Loschmidt echo M n . In the case of S n and M n we measure the approach to equilibrium and for 
tah purpose we explicitly subtracted the invariant state poo in order to avoid saturation. 




4 8 „ 12 16 20 

n 



Figure 3: Time evolution of the linear entropy S n (diamonds) and the autocorrelation function 
(triangles). The points result from the average of ten initial conditions corresponding to coherent 
states centered at ten different points. The Dot-Dash-Dot line indicates the initial Lyapunov 
behavior and the Dashed lines indicate the RP behavior. The slopes of the RP lines were obtained 
from the leading eignvalue computed using the iteration method. As it was expected, the slope for 
Sn is exactly two times the slope of C n . The dotted line shows the behavior if the uniform density 
is not subtracted. 

The long time linear regime can be best observed in Fig. El where C n and S n are plotted. 
In order to uncover the long time regime for S n we subtracted from the initial states. 
The initial Lyapunov regime is very well observed for S n and for longer times the RP regime 
is also well denned. The fact that C n has no Lyapunov regime is well understood. For a 
chaotic system the overlap between the initial state and the succesive iterates is a highly 
fluctuating function until after the Eherenfest time (in which the phase space distribution 
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has reached the edges of phase space). This is the reason why we averaged over a number of 
initial states in order to get a smoother function. The time graphs for the Loschmidt echo 
are very similar 10 



IX. CONCLUDING REMARKS 



Using some new numerical methods we established quantum classical correspondence 
for chaotic quantum maps throught the spectrum of the coarse-grained propagator. For a 
wide range of parameters, the spectrum could be related to the classical Ruelle Pollicott 
resonances. Many techniques to compute the spectrum where presented and some of the 
advantages of the iteration method (mainly for chaotic systems) could be appreaciated. 
As an open question there remains the issue of establishing an optimal relation betwwen 
the noise strength e and h (N) so that the classical limit can be attainedcorrectly. Recent 
work[60] provides an expression such relation in the context of correspondence-breaking time. 
Although it looks promising we have no conclusive proof of its relation to the spectrum. 
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